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Abstract 

Quadrupole radiation of an atom in an arbitrary environment is investigated within classical as 
well as quantum electrodynamical approaches. Analytical expressions for decay rates are obtained 
in terms of Green's function of Maxwell equations. The equivalence of both approaches is shown. 
General expressions are applied to analyze the quadrupole decay rate of an atom placed between 
two half spaces with arbitrary dielectric constant. It is shown that in the case when the atom is 
close to the surface, the total decay rate is inversely proportional to the fifth power of distance 
between an atom and a plane interface. 
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I. INTRODUCTION 



In the recent years the goal of science now has been moving from consideration of the 
fundamental properties of atoms to controlling and changing these properties. It is well 
known that decay rates of atoms can be changed in cavities Many investigations, have 
been devoted to a description of the cavity QED effects Q|. The main attention was paid 
there to the allowed (dipole) transitions. The decay rates of dipole transitions were investi- 
gated in the vicinity of spherical, cylinder, cone, spheroid, aperture, and more complicated 
nanobodies 

However, the influence of environment on the forbidden (quadrupole) transitions is also 
of great interest. First of all, it can help one to study the forbidden transitions. Second, 
with the help of the forbidden transitions one can describe the long-living states, which are, 
in turn, very important in many applications (quantum computers, quantum information). 
Finally, the atoms or molecules with forbidden (quadrupole) transitions can be used as 
detectors of field inhomogeneites. 

The first experiment dealing with quadrupole transitions near a plane interface was car- 
ried out not along ago. The influence of interface on the absorption of Cs 6 2 Si/ 2 — 5 2 D 5 / 2 
transition was studied j^, f^Q. As for the theoretical works, there were very few analyses 
on this topic. In the classical calculations of decay rates of quadrupole transitions 

near the plane dielectric interface were performed. The quadrupole transitions near sphere 
and cylinder were considered within both the classical and QED approaches, and it was 
shown that both approaches gave identical results, as shown in [l(|>[n| . However, there 
was no exact proof of equivalence between the classical and QED pictures. 

The aim of this paper is to find expressions for the quadrupole decay rates and to prove 
their equivalence in an arbitrary environment. In Section |H] we derive expressions for total 
decay rate of a Lorenz oscillator in arbitrary environment through Green function of Maxwell 
equations. In Section ITTT1 we find the expressions for total decay rate of an atom in arbitrary 
environment within the QED approach, and show that they are the same as those in the 
classical approach. Then we apply general results to find expressions for quadrupole decay 
rates for an atom placed in a planar cavity (Section IIV|) This problem is very important 



for some experiments on reflection spectroscopy in thin cells 



12|. General expressions for 



quadrupole decay rates in planar cavity are investigated for the case of dielectric or metallic 
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walls in Section IVl 



II. CLASSICAL DESCRIPTION OF QUADRUPOLE RADIATION IN ANY EN- 
VIRONMENT 



yn a 
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It is easy to show within classical electrodynamics that the total rate of work performed 
by the field can be presented in the form of 



^= / d 3 rJ(r,t)-E(r,t) 



(1) 



where J (r, t) and E (r, t) are the density of current and strength of the electric field, re- 
spectively. This power represents a conversion of electromagnetic field into mechanical or 
thermal energy. In quasi monochromatic case we have instead of the following expression 



dE 



^Re J d 3 vJ* (r, u) ■ E (r, uo) 



(2) 



where u is the frequency, and * means the complex conjugation. 

The electric field can be expressed through current with the help of retarded Green 
function: 



Ei (r; u) = % - J dYGg (r, r'; u) J, (r' ; u) 



(3) 



Here and below the lower Latin subscripts denote Cartesian coordinates and are to be 
summed over when repeated. 

The retarded Green function Q is the solution of Helmholtz wave equation 



V x ( V x *G R (r, r'; w)) - 



e(r)G R (r,r>:u) 



■ { ^ ~T<Sv 



(4) 



where e (r) stands for dielectric constant of environment, and for simplicity we assume that 
the media are nonmagnetic and nondispersive. 

Substituting in (J2J the expression for power can be presented in the form 



IF = ~L Im I dh I d3r ' J * (r ' u) G ^ (r ' r ' ; u) - h (r ' ; u) 



(5) 
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To compare classical and quantum calculations it is convenient to consider stored energy 
as £0 = faj- As a result the expression for decay rate will take the following form 



For the relative decay rate we will have respectively 
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class Im j rf 3 r J d 3 r 'j* fa ^ QR fa J. ^ j. fa. 



(7) 

ittfi Im j dH j dH'J* (r, u) G§° (r, r'; u>) J, (r' ; cu) 
where ^f^o and Gf-° (r, r'; uS) are the total decay rate and the Green function in uniform 
(free) space, respectively, 



Gi- (r,r,w) 



k 2 (5 i:j - nuij) 1 + (3n ini - 6ij) ( - — ^ - - 2 ) 
|r — r | \ r — r r — r / 



gifc|r— r' 



r — r' 



(8) 

In (IHI), n = - r is the unit vector in the direction from the atom to the observation 

|r-r'| 

point and = u/c is the wave vector of free space. 

In the case of quadrupoles the current and charge densities have the following form 



~X)Q<iV i V i *(r-r') (9) 



3! 
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j? = |f£<w(r-ro ( 10 ) 

where r is the radius- vector of the observation point, r' is the radius- vector of the 
quadrupole position, and Qij is the traceless quadrupole momentum tensor 

Qij = J rfxp(x) (3xiXj - x 2 5ij) . (11) 

As is known [l^, any quadrupole can be built out of two dipoles that are equal in 
amplitude and opposite in orientation. In Fig|T]some quadrupoles and the related quadrupole 
momenta are shown. 

Now, after a substitution of (jl())l into (JBJ) and ((Zj) and a partial integration, the final 
expressions for full and relative decay rates will take the following form 

lZT S ' Q = ^ Im ^QhQkNi Vfil (r, r'; W ) (12) 
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dass,Q MKmQJ i Q JH V i V / { G& (r, r'; u) 

I tot r^r' 



(13) 



7t t,o ' Q Im limg*.g w V,V;Gf/ (r, r'; W ) 

r— >r' 

In ()12)) .([13| ) and hereafter, V, V' mean the differentiation over r or r/ , respectively. By 
calculating the limit in ()12j) for free-space Green function (|SJ), we obtain that the expression 
for quadrupole decay rate in free space gets the following simple form 

i»*=3ssEw«i J (14) 

Substituting this expression into (fT3*j) we obtain the following expression for relative 
quadrupole decay rate 

dass,Q Im limQf .QjyVjVjG^ (r, r'; w) 
^f-^=5^^ 5 (15) 

tS,o' q fc 5 El^il 2 

Thus, to calculate quadrupole decay rates it is suffice to determine the respective deriva- 
tives of Green function of Helmholtz wave equation. This result was quite expected as we 
know that electric quadrupoles interact with field gradients. 



III. QED DESCRIPTION OF QUADRUPOLE RADIATION IN ANY ENVIRON- 
MENT: LINEAR RESPONSE THEORY 

To calculate the quadrupole decay rate in arbitrary environment we use the work of 
ref. 3 j 13' but apply a minimal coupling Hamiltonian with generalized Coulomb gauge 

2 

if = — ( pA (r) + Ap (r) ) + -^-A 2 

2mcV y ! PV V 2mc 2 (16) 

div (e (r) A (r)) =0, <p = 

Here p is the operator of electron linear momentum and A (r) is the vector potential at 
the electron position r. The last term in (|16J) gives no contribution to calculation of decay 
rates. 

Assuming that the matrix element of the electron momentum between initial |i)and final 
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\f) states is zero, that is, (/|p|«) = , the Hamiltonian (|16p can be presented in the form 

Hmt = ~^c^ Al (r,) iPl {rj ~ r ' j) + {rj ~ T ' j) Pl) 
where r' is the vector of atom position. 

In first order , the transition rate from initial atomic state |z)to a final state |/) is given 
by Fermi's golden rule 

= t (r) (r ° |/)x (is) 

5 (E F + E } -E I - Ek) 

where = (f\ ((r 3 - — r'j)pi) |i)and capital letters denote eigenstates of the rest of the total 
system under consideration, neglecting its interaction with the atom of interest. Such eigen- 
states might involve, and depend on the coupling between the radiation field, other atoms, 
surface excitations, and the like. For convenience, we refer to those as the "field states". 
For simplicity we assume here that the field is in thermal equilibrium at a temperature T; 

p(I) = exp(— PEj)/Y^ exp(— /3E k ) with j3 = (fc^T) -1 , is the probability that the field is in 

k 

state I. 

Expressing the 5 function of (fTH|) in the integral form we find 

oo 
— oo 

where ujq = (Ef — Ei) /h. In Eq. (119)1 angular brackets indicate an ensemble average and 
A (r ,t) is an interaction picture operator, evolving as if (|16|) were not present, 

A (r, t) = exp (jjH t\ A (r, t) exp ( ~H tj (20) 

In (120)1 is the Hamiltonian of the whole system without interaction. 
Further, one can rewrite (|19j) as Fourier component of the two-point correlation function, 
G? i ,( T y;t) = (A i (r,t)A, (r',0)) 

R " = f (ryi ^ (21) 
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As in our gauge E = — it is possible to show, that ({21} can be presented as 

** -pfe)' (r/; ^ « (22) 

where Gf it (r, r'; t) = (Ei (r, t) (r', 0))is the two-point correlation function of electric field. 
It is convenient to express it through the retarded Green function defined as 



G^(r,r';0 = ^([^(r,*)^(r',0)])e(t) (23) 



In ()23p square brackets mean a commutator and G 
By applying the fluctuation-dissipation theorem 



t) is the Heaviside step function, 
we obtain 



p _ 1 ( e \ ,: m 9 9 nf*n*f 2/ « r - r ^0) (OA) 

Kft - h \mujj SSdr, Or'y >' yv [1 - exp(-f3hw ] 1 j 



where 



oc 

Gfi, (r, r'; cu ) = J dtG? v (r, r'; t) exp (iu t) (25) 

— oo 

is the Fourier component of retarded Green function of electric field. 

The temperature dependence, which appears in the form of an occupation number will be 
important only for (&bT) > Iuoq. Since we are interested primarily in the atomic transition 
energies of the order of a Rydberg, we can set T = OK in this equation. As a result the 
quadrupole decay rate will have the following form for T = 0: 

% = \ (—) ^r^4r DS jl D %i ImG * ( r > r > "o) (26) 

n \mujQ j r^r' or j or-, J J 
As G^, (r, r'; uS) describes the response of the system, it is possible to show that this 
function is the solution of Maxwell equations : 



Vx(vx?(r,r'; W ))-(^ 2 £ (r) 5* (r, r'; v) = 4tt * T<5 (r-r') (27) 



J 2 e(r)C?(r,r»=47r * 

For quadrupole transitions with changing of principal or orbital quantum number the 
following identity is true 



d 



dxj 



^</I^K> (28) 
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Substituting it into eq ({231 and using definition of quadrupole momentum Q{ % - = 
e{(3xiXj — x 2 5ij))j i we obtain the following expression for decay rates for arbitrary 
quadrupole transitions 

R% = ^}™yMQ f Ml™G? k (r, r'; u ) (29) 

where Gf k (r, r'; u Q ) is the retarded Green function of Maxwell equation (|27jl . 

It is very important to remember that this expression is valid for any media, including 
media with losses. 

The quadrupole decay rate in free space is described by the same expression but with 
free space of the Green function (r, r'; u> ), instead 

Rlf = ^}™yMQMk ImG %° ( r > r > «o) = i^'i 2 ( 30 ) 

As a result relative decay rates gets the following form 

R Q MQ^ImGg (r, r'; u ) 

4'° ^ 5 EIQ,l 2 1 ] 

Comparing expression (fM|) with the classical expression (fTo"]) one can see that they are 
identical. It means that both classical and QED models are equivalent for description of 
the total decay rate. Comparison of ()30|) and (fPfj) reveals the difference by the factor of 
4. The same difference takes place in the case of dipole transitions and is related to different 
definitions of dipole and quadrupole momenta in classical and quantum mechanics. 

One should also remember that these equations describe the total decay rates. To find 
the radiative decay rates one should use other approaches, which allow one to take into 
account the radiation patterns of photons. It can be done, for example, within the classical 
approach. 



IV. QUADRUPOLE DECAY RATES OF AN ATOM PLACED IN A PLANAR 
CAVITY 

To calculate the decay rates of quadrupole transition in an atom placed between two 
dielectric half-spaces (Fig. El ) one should find the electric Green function of Maxwell equa- 
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tion. It is very important to ensure that this function should satisfy the symmetry condition 
and the Lorenz reciprocity relation, which follows from the definition (1231) . The approach 

n n 

suggested in [19( allows one to build such a function. According to [19] the Green function 
in layered media can be presented in the following form (z > z') 



d 2 k 



G (r, r'; u) = j ^p-^G (k 



G (k, z, z'; uj) 



(2tt)' 
2m k\ 



E>{k,u,z)E<{ -k,u,z') 



q=p,s 



gl 

1 'I2'l3 e 



(32) 



In fj32|) . (k, u>, z) , i?5 ( — k, u, z') are the mode functions 



E> (k, u, z) = e+ (k) e^(-£) + r ? 2 eC ( k ) e 



E< (k, W| z') 



e^i (k) e~ i/3 ^ + rf 3 ej (k) e i/3 i 2 



and 



(33) 



(k) 



fci 



e± (k) = k x z 



TPik + kz 



- P i 



(-k) 



"si 



-k) 



(34) 



Here = ^Jk 2 — k 2 = y 1 £jk\ — k 2 (k — u/c) is the longitudinal wave vector and r^ 2 , 
are the conventional Fresnel reflection coefficients 



e j P i + e i j ,rij 



Pi " P* 



(35) 



for p and s polarized waves, and L is the distance between plane interfaces. 

Now by substituting this function into ()29j) and integrating it over angle p> in x-y plane 
(k x = cos (f, k y = simp) we obtain the expression of the quadrupole decay rate. 

In the case of the z-oriented quadrupole, that is in the case when 



-1/2 
Q = Qzz -1/2 
1 

the decay rate in free space according to (fHUj) gets the following form 



(36) 



7° 

I zz 



60 



(37) 
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and the expression for relative decay rate has the following form 



To ZZ ~M J 





(38) 



(1 - r\ 2 rl,e^) 

In the case of a single interface with ((L — 2 ) — s — ¥ °°) we have a more simple result 



{£) = W Re I k ' dk(31 (1 " r i3 e24/3l2 °) = 1 - ^ Re J k 3 dkM 3 e 2i ^ (39) 
zz 

This coincidence is very interesting because the Green function used for calculation of 

decay ratej^] is asymmetric. 

In the case of xy+yx quadrupole or in the case of xx-yy quadrupole, where 



Q = Q 



■ni 



1 

1 




1 

Q = Qxx o-io 
ooo 

the decay rates in free space according to (jSOjl get the following form 



(40) 



(41) 



7 
7 

and for the relative decay rate we have, respectively, 



_ 




xy 




_ 




XX 


45 



\ Q 

-) 

7o/ xy 



J_ 

7o / a ; 



oo 

5 f k 3 dk 

M ' J ~K 





2 (l-^ 2 e^^)(l-^ 3e ^i-0) 
Pi (i_ r P P 2H3 lL ) 



(42) 



(43) 



(l-rJ 2 rf 3 e : W) 

In the case of s — > oo, that is in the case of single interface, we have a more simple result 
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= 4^/ ^ [Pi (1 - r p 13 e^) +*?(! + rfee^*)] 



o ft 
°°k 3 dk 



2i/3 1 za 



Finally, in the case of xz+zx or yz+zy quadrupoles, where 



(44) 



Q Qxz 



Q = Q 



1 



1 



1 

1 



the decay rates in free space according to (J30|) get the following form 



k^Q 



2 

xz 



k()Qyz 



ijz 



45 



and for relative decay rate we respectively have 



(45) 



(46) 



(47) 



j_ 

7q/ a : 



1_ 

70/ y: 



oo 

5 f kdk 

Ml J "ft" 





, 2 _ 2 x2 (l+rg 2e ^i-)(l+rf 3 e^i'0) 

{Pi K ) (l-^ 2 rf 3 e 2 ^i L ) 

n2, 2 ( 1 - r ^ 2il3ls )^-rhe 2i ^o) 
Pl K l d_ r s s 2 t ^ Ll \ 



+ 



(48) 



In the case s — > oo, that is, in the case of single interface, we have a more simple result 



kdk 



4fcn 5 "~ J ft L 





5 „ ™kdk r,„, ,on2 



1 + 4F^1T M-* 2 ) ^-'"^ 

4K Pi 



O 2i/3 1 z 



(49) 



As mentioned above, these results describe the total decay rates, i.e. radiative and 
nonradiative. Generally, it is difficult to separate these contributions. However, one can 
assume that this separation can be made on the basis of the classical energy flux method. 
It should be noted that a purely radiation channel may exist, in this geometry, in an ideal 
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case of matter without losses. At negligibly small losses, the radiation energy would not go 
to infinity. This is the difference between the geometry under consideration and an open 
geometry, at which the radiation might go to infinity throughout a free space. 



V. ANALYSIS OF RESULTS AND ILLUSTRATIONS 

The expressions that had been obtained in the previous section are rather complicated 
and their calculation is an independent problem, in a general case. The complexity is due to 
the fact that the integrands are the complex functions with a set of the singular points, which 
might be both the branching points, and the poles. These peculiarities are connected with 
physical properties of the problem. In any case, in the integrand there are the branching 
points at k — ±y/eko- If the mode wave propagation is formed in a cavity (metallic mirrors) 
then the poles appear in the integrands. So, in different physical situations, the calculations 
are to be performed with account of these factors. 

A. Atom between perfect metallic mirrors 

In the case of the well conducting metallic mirrors, the expressions (|38|). (H2J), and 
in which the reflection coefficients are substituted by their analogs for the case of an ideal 
conductivity, 



will be good approximations for the rates 

The expressions for the decay rates of quadrupole states may be reduced to the form 



l,r s 



-1 



(50) 




(51) 



o 



Q 



Q 



oo 




.''.V 



7o 



7 



) 



XX 




sin (Pis) sin (/J^o) 
sin (PiL) 



(52) 



1 ' 
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1_ 

v70/ a-: 



_7_ 

70/ y: 



2&0 5 ""7 £1 





cos (/^s) COS (P1Z0) 
sin (/^L) 



(53) 



Because the ideal conductivity is the limiting case for a real metal, where the poles must 
lie above the horizontal axis of integration, the integration circuit of the ideal conductivity 
must envelope the poles from below, as shown in FigOJl 



To calculate the integrals ()51| ) -([53 J) it is convenient to use the variable /3 X = \fk$ — k 2 , 



Q k o 
7o / zz K o J 



7o / xy 



feo 



_7 

70/ a-.r 



sin (/3 1 s) sin (fliZo) 
sin(^L) 

sin (/?is) sin (fliZo) 



(54) 



(55) 



-) 



_7 

70/ j/2 



Im / dp x (201 - klf + 0{kl 



COS (/^s) COS (@iZq) 

sin (/^L) 



(56) 



where the path of integration is shown in FigJS]. By calculating the integrals (|54j ) -([56j ) with 
the residue theorem one can obtain the following results 



j_ 

70/ zz 

-) Q - 

70/ xy 



1 r- ^max / \ 2 

f £ (?) ( 1 



7m 



57T 



n=l 



I 



. 2 /7mz \ 



8111 ( — 



(57) 



(58) 



\ Q 

-) 

7o/ 



where n r 



\ Q 

-) 

7o/ yz 



5tt 
2L 



-1 %ax I 

5 + E 



n=l 



2, T 



7m 
T 



2 /7TnZ \ 



(59) 



= L/7T is integral part of L/ir, and L = k L 
Figure 0] illustrates the quadrupole decay rates in a resonator formed by a hypothetic 
metal with e = — 200+0. Oli in respect to the position and orientation of a quadrupole. 
As seen from the Figure, the asymptotic expressions (f5T J) -([5l? |) approximate well the exact 
expressions (J3*8|) . (J4*3~)) (}4*Sj) . excluding the region that is in a close proximity to the metal 
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surface. However, in the vicinity of the surface, the nonradiative losses connected with 
imaginary part of the dielectric constant are of the main importance. These losses are the 
reason of a fast increase in the total losses (see Eqs.([64 |) -([66 )l ). In the case of an ideal 
conductor, the losses are absent, and there is a difference between the decay rates of a 
hypothetic metal and an ideal conductor. 

In the case of real metals, that difference might be still more profound because the 
imaginary part of the permittivity is not negligibly small as compared to the real part. 
Figure 03 illustrates the decay rates for a micro-resonator with silver mirrors. From the 
Figure one can see that the rate of spontaneous decays in the real resonator differs from the 
decay rate in the cavity with ideal walls substantially. 

B. Atom between dielectric mirrors 

A planar cavity can also be realized on the basis of two opposite dielectric half-spaces. No 
propagating waveguiding modes are formed in that case, and the integrand, respectively, has 
no poles in a complex plane near a real axis. This should simplify a numerical calculation 
of the integrals. Figures 16171 demonstrate the dependencies of the quadrupole decay rates 
on the atomic position and structure of the quadrupole moment for a planar resonator with 
silica walls. 

In the case of quite a large-size cavity (micro-cavity, FigEJ) one can observe an increase 
in the rate of spontaneous decays as an atom is approaching the wall. In contrast to the 
case of metallic mirrors, such an increase is due to the coupling of the non-propagating near 
fields emitted by the quadrupole, with the propagating fields inside the dielectric (silica). 
One can notice the influence of the intrinsic non-radiative processes at a distance less than 
1 nm only, because the imaginary part of the quartz permittivity is very small at optical 
frequencies (see Eqs. dOIJ)-(JEEJ)). At such distances one should take into account the random 
inhomogeneities of the surface structure. 

As the distance between dielectric walls is small (nano-cavity, Fig|7j), the electric fields 
are near fields one, at any atom position between the walls, and there occurs the effective 
field transformation into the wave propagation over a dielectric. This provides a consider- 
able acceleration of the transitions. The intrinsic non-radiative decay channel is formed at 
distances closer to the wall, and this is unseen on the picture. 
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Note that all the energy of an excited atom will be emitted in the dielectric, and all the 
losses will, therefore, be non-radiative. But in the case of the weakly absorbing dielectrics, 
including silica, it is not unreasonable to distinguish between the regions of the effective 
transformation into the propagating waves and the regions of the intrinsic radiative losses. 



C. Atom inside ultra thin cell 



Very interesting spectroscopy experiments are carried out now with atoms inside an extra- 
thin dielectric cell. Suffice it to say that the width of the cell cab be as small as 20 nm . 
So, it is very interesting to understand the behavior of decay rate in that case. 

All dimensional parameters are small in comparison with wavelength. As a result we 
can use the quasi-static approximation to calculate decay rates found in a previous section. 
The quasistatic approach here is equivalent to the case of k >> k = uo/c. In this limit the 
Fresnel reflection coefficients can be simplified substantially 



p v 

r 12 = r 13= r 



la = r t3 = 



E-l 
6 + 1 



(60) 



As a result the decay rate in the small width of cavity case will have the following form: 



j_ 

70/ zz 

\ Q 

-) 

7o/ xy 



45 x 
Im > 



( 



2n+l 



n=0 



2r 



1 



L(n + 1) 



15 
16 



/ 



Im J2r 2n+1 



n=0 



2r 



/ 



v 



L(n + 1) 
2r 



L(n+1) 



+ 



s + Ln 



s + Ln 



s + Ln 



+ 



zq + Ln 



1 



Zq + Ln 



Zq + Ln 



(61) 



(62) 



(63) 



/ 



where z ,s,L stand for k z , k s, k L , respectively. In the case when atom is very close to 
one surface Zq << L only one term ( n=0 ) is important in this series 



j_ 

7o 



Q 



45 

si5 



Imr 



(64) 



15 





^ 15 


Toy 






\ ® 15 


Toy 





Imr (65) 



Jmr (66) 

From this asymptotics one can see that the total decay rate increases inversely propor- 
tional to the fifth power of distance to surface z — 0. This behavior is different substantially 
from the dipole case, where decay rates increase inversely proportional to the third power 
of distance to surface. Another interesting point one can get from (|64j) - (JBTIj) , is that the 
zz-quadrupoles suffer a six-fold enhancement in comparison with other components. 



VI. CONCLUSIONS 



In this article, the processes of the spontaneous quadrupole atomic radiation in an ar- 
bitrary environment were considered within the framework of both classical and quantum 
electrodynamics. The general equations derived for the rates of quadrupole transitions were 
expressed through the spatial derivatives of the retarded Green function corresponding to 
the classical problem of electrodynamics. It was shown that the expressions differ by a 
numerical coefficient 4 only, which is connected with the definition of quadrupole moments 
which have different physical sense in the classical and quantum mechanics. The expressions 
for the relative decay rates, i.e. the rates normalized by the uniform space rate, prove to be 
identical. 

The results obtained are applied to a description of quadrupole atomic transitions in a 
planar cavity. The explicit analytical expressions for the rates of any quadrupole transition 
were found for such a cavity. The results have been analyzed in detail for the planar 
cavities with dielectric and metallic walls. It was found that the quadrupole transitions are 
accelerated with decreasing resonator size. In the case of dielectric walls, such an acceleration 
is due to the transformation of the near dipole fields into the propagating waves inside the 
dielectric. In the case of metallic mirrors, the acceleration becomes more profound, and is 
due to the radiation absorption at the surface layer of a metal. 

In this paper we restrict ourselves to investigation of quadrupole decay rates. However, 
our approach can be also applied to description of frequency shifts of quadrupole transi- 
tions in nanoenviroment. Again, general expressions for frequency shift will be expressed 
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through space derivatives of retarded Green function. We will present detailed investigation 
of frequency shifts of quadrupole transitions in nanoenviroment in a separate publication. 
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FIG. 1: Some kinds of quadrupoles and the corresponding quadrupole momentum tensors. 



19 



z 









k A 


S 

^Quadrupole 
Z 8, = 1 

r 







FIG. 2: Geometry of the problem of quadrupole radiation of an atom placed in a planar cavity. 
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FIG. 3: Contours of integration in the case of ideally conducting mirrors. 
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FIG. 4: The quadrupole decay rates of different quadrupoles versus their position in the case of a 
hypothetic material with e = —200 + O.li, and in the case of an ideal conductivity (dotted line). 
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FIG. 5: The quadrupole decay rates of different quadrupoles versus their position between two 



thick silver (Ag : e = 15.37 + 0.231* ,A 
case of the ideally conducting walls). 
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FIG. 6: The spontaneous decay rates of different quadrupoles versus their position between two 
quartz half-spaces with e = 2.1 + O.OOOOOOOOli(silica) in the case of a micro-resonator (kL = 10). 
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FIG. 7: The spontaneous decay rates of different quadrupoles versus their position between two 
quartz half-spaces with e = 2.1 + O.OOOOOOOOli(silica) in the case of a nano-resonator ( kL = 1). 
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